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is discussed and the initial values of the doubled variables are related to quantum noise 
effects. 
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1. Introduction 



A dissipative system is physically incomplete and a microscopic theory must 
include the details of processes responsible for dissipation, including quantum ef- 
fects. One then would start from the beginning with a Hamiltonian that describes 
a complete system, the bath and the system-bath interaction. Subsequently, the 
description of the original dissipative system is recovered by the reduced density 
matrix obtained by eliminating the bath variables which originate the damping 
and the fluctuations. The problem with dissipative systems in quantum mechan- 
ics is indeed that canonical commutation relations (ccr) are not preserved by time 
evolution due to damping terms. The role of fluctuating forces is in fact the one 
of preserving the canonical structure. 

The above strategy, however, may not always be viable since it requires the 
knowledge of the details of the processes inducing the dissipation; these details may 
not be explicitely known and the dissipation mechanisms are sometime globally 
described by such parameters as friction, resistance, viscosity etc.. In such a case a 
different route must be followed. On the other hand, the attempt to derive, from 
a variational principle, the equations of motion defining the dissipative system 
requires the introduction of additional complementary equationsfl]. 

The latter approach has been pursued in refs. 2 and 3 where the quantization 
of the one-dimesional damped linear harmonic oscillator (dho) has been studied 
by doubling the phase-space degrees of freedom. The new degrees of freedom thus 
introduced play the role of the bath degrees of freedom and make it possible to 
derive the dho equation from a variational principle. 

The quantum hamiltonian is obtained as 



where T is the decay constant, O the frequency and A\ , A and B are creation 
and annihilation operators with usual ccr [2, 3]. 

As pointed out in refs. 2 and 3, only if one restricts the Hilbert space to 
the subspace of those states, \ip >, which are annihilated by B: B\ip >= 0, the 
eigenstates of H merge into those of the simple undamped harmonic oscillator 
when r — > 0. Therefore the states generated by B^ represent the sink where the 
energy dissipated by the quantum damped oscillator flows: the S-oscillator thus 
represents the reservoir or heat bath coupled to the A-oscillator. 

The dynamical group structure is that of SU(1, 1) and it has also been 
shown [3] that the canonical quantization of the dho can be achieved in a consistent 



H Q = hQ(A*A - B^B) 



H = Ti + Tii i 
, Hi = iUY(A^B^ - AB) 



(1) 
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way in the infinite volume (or thermodynamic) limit of Quantum Mechanics, i.e. 
in Quantum Field Theory (QFT) where infinitely many unitarily inequivalent rep- 
resentations of the ccr are allowed. The reason for this is that the set of states of 
the dho splits into unitarily inequivalent representations (i.e. into disjoint folia, in 
the C* -algebra formalism) each one representing the states of the oscillator at time 
t and therefore QFT provides the proper setting to discuss dho: the irreversible, 
non-unitary time evolution is thus described as tunneling across unitarily inequiv- 
alent representations (breakdown of time reversal invariance or arrow of time) . In 
such a description a central role is played by the degeneracy among the vacua of 
the inequivalent representations, which is formally expressed by the fact that the 
free Hamiltonian H is the Casimir operator of the SU{1, 1) group and therefore 
a constant of motion, [7i ,Ti.i] = 0. Quantum processes changing separately ua 
and n b are allowed provided their difference is left constant and therefore H may 
have (infinitely) many degenerate zero-energy states (vacua). The positiveness of 
H is also ensured by its constancy in time, once its initial (positive) value has 
been assigned. 

In ref. 3 it has been shown that the dho states are time dependent thermal 
states, as expected due to the statistical nature of dissipation. In particular the 
formalism for the dho turns out to be similar to the one of real time QFT at finite 
temperature, also called thermo-field dynamics (TFD)[4], which is indeed based 
on the doubling of the degrees of freedom. 

The purpose of this paper is the following. In the approach of refs.1-3 the 
system is perturbed in a manner that its time development is governed by different 
dynamics depending on the positive and negative sense of time. In this work we 
want to show that there is consistency between the variational principle in the 
Lagrangian formalism and such different dynamical behaviors relative to a different 
sense of time. Or, in other words, on a closed time path the variational principle is 
sufficient to determine the expectation values of any physical observable for a given 
initial condition. We will obtain the exact action for the dho in the path integral 
formalism. In particular we will reach some conclusion on the initial values of 
the doubled variables and relate them to the probability of quantum fluctuations 
in the ground state, a result which is interesting also in the more general case of 
thermal field theories, such as TFD, with the doubled variable formalism. 

We will cast the canonical quantization scheme of the dho of ref. 3 in the 
formalism developed by Schwinger[5] and Feynman and Vernon[6] since it is par- 
ticularly suited to our task. 

2. Lagrangian formalism 
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Our purpose is to explore the manner in which the Lagrangian model for quan- 
tum dissipation of refs. 2-3 arises from the formulation of the quantum Brownian 
motion problem as described by Schwinger[5] and by Feynman and Vernon [6]. Let 
us first recall (and generalize) the Lagrangian of dho[2,3]. The Langrangian for a 
particle of mass p,, damped by a mechanical resistance R moving in a potential V, 
reads 

C(x, y, x, y) = pxy - V(x + (l/2)y) + V(x - (l/2)y) + (R/2) (xy - yx). (2) 
The classical Lagrangian equations of motion implied by Eq.(2) read 

fix + Rx+ (1/2) [V'(x + (l/2)y) + V'(x- (l/2)y)] = 0, (3a) 

fiy-Ry+ [V'(x + (l/2)y) - V'(x - (l/2)y)} = 0. (36) 

It is easy to see that for V(x ± (l/2)y) = (l/2)k(x ± (l/2)y) 2 Eqs. (3) give 
the dho equation and its complementary equation of refs. 2 and 3. 

If from the manifold of solutions to Eqs. (3) we choose those for which the y 
coordinate is constrained to be zero, then Eqs. (3) simplify to 

y = 0; fix + Rx + V'(x) = 0. (3c) 

Thus we obtain a classical damped equation of motion from a Lagrangian theory 
at the expense of introducing an "extra" coordinate y, later constrained to vanish. 
(Note that y(t) = is a true solution to Eqs. (2) so that the constraint is not in 
violation of the equations of motion.) 

The question to which we address ourselves is the following. Does the intro- 
duction of an "extra coordinate" make any sense in the context of conventional 
quantum mechanics? Consider the special case of zero mechanical resistance. For 
that case, one should begin with the Hamiltonian for an isolated particle, 

H = -(h 2 /2/j)(d/dQ) 2 + V(Q). (4) 

For an isolated particle one obtains in quantum mechanical theory the density 
matrix equation of motion 

ih(dp/dt) = [H,p], (5) 

which indeed requires two coordinates (say Q+ and Q-). In the coordinate repre- 
sentation, Eqs. (4) and (5) read 

ih(d/dt) < Q+\p(t)\Q- >= 
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{-(h 2 /2p)[(d/dQ + ) 2 - (d/dQ_) 2 ] + [V(Q+) - V(Q_)}} < Q+\p(t)\Q. > . (6) 
Employing the coordinates x and y, 



Q± = x±(l/2)y, (7) 

and the associated density matrix function 

W(x, y,t) =<x + (l/2)y\p(t)\x - (l/2)y >, (8) 

simplifies Eq.(6) which now reads 

ih(d/dt)W = H W (9a) 

n o = ( PxPy /fi) + V(x + (l/2)y) - V(x - (l/2)y), (96) 
p x = -ih(d/dx), p y = -ih(d/dy). (9c) 
Of course the "Hamiltonian" Eq.(9b) may be constructed from the "Lagrangian" 

C = fixy - V(x + (1/2)2/) + V(x - (l/2)y). (10) 

We have then the justification for introducing Eq.(2) at least for the case R = 0. 

Finally we observe that for V(x ± (l/2)y) = (l/2)k(x ± (l/2)y) 2 it is easy to 
obtain the Hamiltonian (1) from the Lagrangian Eq.(2)[2,3]. 

We also notice that 7i and Tij in Eq.(l) are the free Hamiltonian and the 
generator of Bogolubov transformations, respectively, in TFD[4]. Our present dis- 
cussion thus includes the doubling of degrees of freedom in finite temperature 
QFT. 



3. Damping 

Now let us suppose that the particle interacts with a thermal bath at temper- 
ature T. The interaction Hamiltonian between the bath and the particle is taken 

as 

Him = ~fQ, (11) 

where Q is the particle coordinate and / is the random force on the particle due 
to the bath. 
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In the Feynman- Vernon [6] quantum Brownian motion scheme, the effective 
action for the particle has the form 



it 

where C Q is defined in Eq.(lO) and 



/•*/ 

A[x,y]= dtC (x,y,x,y)+I[x,y], (12) 



e (i/h)X[ X ,y] =< (WVSl; fWQ-(t)dt)^ e (i/h)f*f f(t)Q + (t)dt) )+ > (l3) 

In Eq.(13): (i)The average is with respect to the thermal bath; (ii) "(.)+" denotes 
time ordering and "(•)-" denotes anti-time ordering; (iii) The c-number coordi- 
nates Q± are defined as in Eq.(7). We observe that if the interaction between the 
bath and the coordinate Q (i.e Hint = —fQ ) were turned off, then the operator 
/ of the bath would develop in time according to f(t) = e lH R t / h j e -iH R t/h w ] iere 
Hr is the Hamiltonian of the isolated bath (decoupled from the coordinate Q). 
/(£) is the force operator of the bath to be used in Eq.(13). 

Assuming that the particle first makes contact with the bath at the initial 
time U, the reduced density matrix function in Eq.(8) obeys at a final time 

/oo poo 
dxi / dyiK(xf,yf,tf;xi,yi,ti)W(xi,yi,ti), (14) 
-oo J — oo 

with the path integral representation for the time development 

rx(tf)=Xf Mtf)=Vf ....... . 

K(x f ,y f ,t f ;x i ,y i ,t i )= / Vx(t) / Vy(t)e^ A ^. (15) 

Jx(ti)=Xi Jy(t l )=y l 

The evaluation of T[x, y] for a linear passive damping thermal bath involves 
several Greens functions all of which have been discussed by Schwinger[5] in his 
work on quantum Brownian motion. For completeness of presentation we here 
discuss these results. 



4. Greens Functions 



The fundamental correlation function for the random force on the particle 
due to the thermal bath is given by 



G(t-s) = (i/h)<f(t)f(s)> 
6 



(16) 



From this one may construct the time ordered Greens function and the anti-time 
ordered Greens function 

G+(t -s) = 9(t - s)G(t -s) + 6(s - t)G(s - t), (17a) 

G_(t -s) = 9(s - t)G(t -s) + 6(t - s)G(s - t), (176) 

where 

9(t- s) = 1 if t > s, (18a) 

and 

$(t- s) = if t < s. (186) 
The retarded and advanced Greens functions are defined by 

G ret (t -s) = 6{t - s)[G(t -s)- G(s - t)], (19a) 

and 

G adv (t -s) = 6(8 - t)[G(s -t)- G(t - s)} . (196) 

From a causal engineering view point, the mechanical impedance Z(Q (ana- 
lytic in the upper half complex frequency plane Xm ( > 0) is determined by the 
retarded Greens function 

/•OO 

-<Z(C)= / dtG ret (t)e iCt . (20) 
Jo 

The time domain quantum noise in the fluctuating random force 

N(t -s) = (1/2) < f(t)f(s) + f(s)f(t) >, (21) 
is distributed in the frequency domain 

N(t-s)= I duS f (u)cos[u(t- s)], (22) 
Jo 

in accordance with the Nyquist theorem 

S f (u) = (hu/n)coth(hw/2kT)neZ(u + i0 + ). (23) 
The mechanical resistance is defined 

R = lim u ^ KeZ(u; +i0 + ). (24) 
7 



Finally, the time ordered and anti-time ordered Greens functions describe 
both the retarded and advanced Greens functions as well as the quantum noise, 

G±(t -s) = ±(l/2)[G re t(t -s) + G adv {t - s)] + (i/h)N(t - s). (25) 



5. Bath Interactions 

For the particle interacting with the bath describing a linear passive mechan- 
ical impedance, Eq.(13) may be evaluated following Feynman and Vernon as, 

I[x, y] = (1/2) 1' T dtds[G+(t - s)Q + (t)Q+(s) 
Jti Jti 

+G_(t - s)Q-(t)Q-(s) - 2G(t - s)Q_(t)Q + (s)} . (26) 
Eq.(26) is more simple when expressed in terms of coordinates x and y, 

l[x,y) = (1/2) 1' 1' dtds[G ret (t - s) + G adv {t - s)][x{t)y(s) + x(s)y{t)} 
Jti Jti 

+(i/2h) f f [ f dtdsN(t - s)y(t)y(s). (27) 
Jti Jti 

Further simplification results if one defines the retarded force on y 

F r y et {t) = J*' dsG ret (t - s)y(s), (28a) 

and the advanced force on x, 

Ff v {t) = I'' dsG adv (t - s)x(s). (286) 
Jti 

The interaction between the bath and the particle is then 

l[x,y] = (1/2) dt[x(t)F^(t) + y(t)F^(t)] 
Jti 

+(i/2h) / / dtdsN(t- s)y(t)y(s). (29) 
Ju Jti 



6. The Final Action 



Putting all the pieces together, we have from Eqs.(12) and (29) the real part 
of the action 



lZeA[x, y] = 




(30a) 



C = »xy - [V(x + (l/2)y) - V(x - (l/2)y)] + (l/2)[xFf* + yF? d % (306) 
and the imaginary part of the action 

TmA[x, y] = (l/2h) / / dtdsN(t - s)y(t)y(s). (30c) 
Jti Ju 

The central Eqs.(30) are rigorously exact for linear passive damping due to the 
bath when the path integral Eq.(15) is employed for the time development of the 
density matrix. 

The lagrangian Eq.(2) can now be viewed as the approximation to Eq.(30b) 
with Fy Ct = Ry and F^ dv = —Rx. The classical constraint y = occurs because 
nonzero y yields an "unlikely process" in view of the large imaginary part of 
the action (in the classical u h — > 0" limit) implicit in Eq.(30c). On the contrary, 
at quantum level nonzero y may allow quantum noise effects arising from the 
imaginary part of the action. 

We are glad to aknowledge useful discussions with E. Celeghini and M.Rasetti. 
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